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ABSTRACT 


-> 


The  present  study  Is  a  benchmark  test  to  aid  the  development  of  various 
rotor  performance  codes.  The  study  Involves  simultaneous  blade  pressure  measure¬ 
ments  and  tip  vortex  surveys.  Measurements  were  made  for  a  wide  range  of  tip 
Mach  numbers  Including  the  transonic  flow  regime.  The  measured  tip  vortex 
strength  and  geometry  permit  effective  blade  loading  predictions  when  used  as 
Input  to  a  prescribed  wake  lifting  surface  code.  It  Is  also  shown  that  with 
proper  Inflow  and  boundary  Tayer  modeling,  the  supercritical  flow  regime  can 
be  accurately  predicted.  ^ 


LIST  OF  SYMBOLS 

A  *  ratio  of  vortex  circulation  to  maximum  blade-bound  circulation 
AR  =  aspect  ratio 

=  sectional  lift  coefficient 

d  *  radial  distance  from  a  vortex  to  a  flow-field  point 
R  =  radius  of  the  rotor  blade 

r  =  radial  distance  from  the  rotor  center  of  rotation 
=  vortex-induced  velocity 

Vjj  =  residual  velocity  in  the  wake 
y  =  r/R,  nondlmenslonal  radial  coordinate 

z  =  axial  distance  from  rotor 

(2  =  rotational  speed 

¥  =  azimuthal  angle  measured  from  the  point  of  blade  overhead  passage 

=  vortex  age,  the  azimuth  angle,  ¥,  when  vortex  strikes  the  probe 

1 .  Introduction 


The  past  two  decades  have  seen  a  continuing  development  of  methods  to 
predict  rotor  hover  performance  with  Inctjeaslng  accuracy.  These  methods 
Include  lifting  line  [1,2,3],  lifting  surface  [4,5,6],  and  finite  difference  [7] 
methods.  Practically  speaking,  none  of  these  methods  Is  self-contained;  they 
all  require  the  specification  of  empirically  obtained  wake  data  (strength  and 
geometry)  In  order  to  have  a  correct  downwash  distribution.  Inevitably,  the 
development  of  these  codes  becomes  a  tuning  process  In  which  It  Is  determined 
Just  how  detailed  and  accurate  a  wake  description  must  be.  This  stage  of  code 
development  places  great  reliance  on  the  available  body  of  experimental  rotor 
data. 


The  available  rotor  data  Include  a  sizeable  number  of  tests  where  detailed 
blade  loading  Is  obtained  using  surface  pressure  transducers  [8,9,10,11],  and 
more  recently  by  laser  doppler  veloclmetry  [12,13].  There  Is  also  a  number  of 
tests  In  which  the  rotor  wake  geometry  Is  defined  by  flow  visualization  tech¬ 
niques  [5,14]  for  a  wide  variety  of  blade  configurations.  Of  the  various  wake 
studies,  only  Boatwright  [15]  and  Cook  [16]  made  detailed  Investigations  of 
the  wake  flow  structures.  Cook's  work  Is  especially  significant  In  that  he 
was  able  to  measure  the  strength  of  the  tip  vortex  by  a  curve-fitting  technique 
using  hot-wire  data.  However,  there  seem  to  be  no  useable  data  In  the  literature 
In  which  simultaneous  blade  load  distribution  and  wake  measurements  are  made. 

It  Is  the  Intention  of  the  present  study  to  help  fill  this  gap  In  the 
literature.  This  paper  will  describe  the  experimental  set-up  In  which  steady 
blade  pressures  were  obtained  using  hub-mounted  transducers  and  tip  vortices 
were  measured  using  Cook's  technique.  The  data  obtained  are  for  unstalled  flow 
ranging  from  the  low  subsonic  to  transonic  conditions.  It  Is  shown  herein  (using 
Summa's  lifting  surface  code  [6])  that  the  blade-loading  distribution  could  not 
have  been  predicted  using  only  the  classical  vortex  data  of  Landgrebe  [14]  and 
Kocurek  [5] . 

2.  The  Experiment 

The  data  presented  In  this  paper  were  gathered  In  the  Army  Aeromechanics 
Laboratory's  hover  test  facility,  a  large  chamber  with  special  ducting  designed 
to  eliminate  room  recirculation.  The  rotor,  situated  In  the  center  of  the 
chamber,  was  mounted  on  a  tall  column  containing  the  drive  shaft  (Fig.  1).  The 
rotor  employed  two  cantilever-mounted,  manually  adjustable  blades.  These  blades 
used  an  NACA  0012  profile  and  were  untwisted  and  untapered.  An  aspect  ratio 
of  6  was  chosen  In  order  to  maximize  Reynolds  Number  and  available  Instrumenta¬ 
tion  space.  The  blades  were  grooved  to  accommodate  60  pressure  tubes  each. 

These  tubes  connect  to  a  special  cluster  of  three  48S8  Scanlvalves  (using 
Statham  PA  856-15  transducers)  driven  by  one  SS5-48  solenoid  drive  mounted  In 
the  rotor  hub.  This  arrangement  permits  an  ample  number  of  ports  for  five 
measurement  locations  —  three  radial  locations  on  each  blade,  with  one  location 
being  Identical  on  both  blades  for  comparison  purposes.  The  Scanlvalve 
stepper  motor  was  actuated  by  a  digital  data  system  which  acquired  the  data, 
computed  the  centrifugal  pressure  drops,  and  displayed  the  final  pressure 
distribution.  The  resulting  pressure  dlstrlubtlons  for  collective  pitch 
settings  of  5°,  8°,  and  12**  are  shown  In  Figs.  2,  3,  and  4.  It  Is  seen  here 
that  the  Inboard  pressure  distributions  are  only  slightly  affected  by  rotor 
speed.  However,  the  outboard  sections  show  considerable  pressure  alteration 
and  shock  development  as  the  tip  Mach  number  approaches  near  sonic  values. 
Overall,  however,  the  spanwlse  load  distribution  (obtained  by  pressure  integra¬ 
tion)  is  remarkably  little  affected  by  tip  Mach  number  (Fig.  5) . 

Wake  data  were  acquired  with  a  traverse-mounted  DISA  55P01  hot-wire  probe 
mounted  beneath  the  rotor.  Data  from  the  wire  are  acquired  at  various  points 
along  the  tip  vortex  trajectories  and  can  give  both  the  tip  vortex  geometry 
and  strength.  One  problem  with  this  approach  Is  that  the  vortex  trajectory  Is 
not  steady  and  the  probe  location  (which  Is  chosen  by  an  on-the-spot  decision 
as  to  where  the  number  of  vortex  core  "hits"  Is  maximized)  contains  some  as  yet 
undetermined  error.  The  resulting  data  stream  has  considerable  variability. 
However,  In  order  to  be  certain  of  the  vortex  location,  the  only  acceptable 
data  are  those  where  the  vortex  core  actually  hits  the  probe.  In  the  digitiza¬ 
tion  process  (done  off-line  at  a  reduced  tape  speed) ,  the  above-mentioned  data 
system  was  coded  to  look  for  and  accept  only  those  data  which  showed  the 
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characteristic  signal  dip  which  Indicates  a  probe-vortex  strike.  This  turns 
out  to  be  a  very  small  percentage  of  the  total  amount  of  data  actually  recorded. 

A  typical  hot-wire  trace  displaying  the  above-mentioned  variability  is  shown  in 
Fig.  6. 

3.  Hot  Wire  Data  Analysis 

The  idea  of  the  current  data  analysis  is  that  a  tip  vortex  should  look 
like  an  infinite  line  vortex  to  a  sufficiently  close  probe.  Unfortunately,  the 
probe  measures  not  only  the  velocity  induced  by  the  vortex  at  hand,  but  also 
that  induced  by  the  blade  and  the  remainder  of  the  wake  system  as  well,  V^. 

The  problem  in  analyzing  the  probe  data  is,  then,  how  to  separate  this  residual 
velocity,  Vp,  from  the  immediate  vortex- induced  velocity,  Vi.  Cook  [16]  handled 
this  probem  by  assuming  that  the  residual  velocity  was  constant  and  given  by  the 
translation  velocity  of  the  tip  vortex.  He  then  was  able  to  find  the  vortex 
strength  by  a  fitting  process.  This  strength  was  found  to  be  far  less  than  the 
computed  maximum  blade  bound  circulation  of  the  single,  full-scale  blade  used  in 
that  test.  It  was  also  found  that  the  vortices  measured  were  distinctly  non- 
classical  in  that  they  contained  a  large  rotational  region  outside  of  the  viscous 
core.  In  what  follows,  we  shall  use  a  process  very  similar  to  Cook's  in  analyz¬ 
ing  wake  data. 

First  consideration  is  given  to  the  vortex  trajectories.  Figure  7  shows 
the  axial  and  radial  components  of  the  vortex  trajectories  for  a  pitch  setting 
of  8°.  This  figure  gives  data  for  a  wide  range  of  rotor  speeds,  and  it  is 
apparent  that  the  trajectory  is  essentially  independent  of  tip  speed  —  even  into 
the  transonic  regime.  Figure  7  together  with  Fig.  5  suggests  that  the  nonlinear 
transonic  flow  on  the  blade  has  little  effect  on  the  far-fleld  induced  flow  as 
long  as  the  local  lift  is  not  greatly  altered.  Also  plotted  on  this  figure  is 
the  vortex  trajectory  given  by  Kocurek's  wake-fitting  formula  for  rotors  in 
free  air.  Although  the  axial  component  of  the  trajectory  compares  well  with 
Kocurek's  formula,  there  appears  to  be  a  greater  discrepancy  in  the  contraction 
than  can  be  explained  by  measurement  error.  The  vortex  trajectories  for  pitch 
settings  ranging  from  5"  to  12°  are  given  in  Fig.  8. 

The  present  aim  in  analyzing  the  rotor  wake  is  only  to  find  the  vortex 
strength  and  not  a  complete  description  of  the  structure.  This  strength  will 
be  found  by  fitting  the  wake  data  to  the  velocities  obtained  from  an  appropriate 
combination  of  invlscid,  two-dimensional  vortices.  The  velocity  from  one  such 
vortex  is  given  by 


^i  ■  <V>max  1 
fiR  “  4ttAR  d/R 


(1) 


where  the  strength  of  the  vortex  is  described  by  A,  the  ratio  of  the  vortex 
circulation  to  the  maximum  bound  circulation  of  the  blade  (determined  from  the 
pressure  data).  To  accomplish  this  fitting,  it  is  first  necessary  to  convert 
the  spatially  dependent  Equation  (1),  into  a  time-dependent  expression,  as  the 
vortex  data  are  time-based.  Assuming  that  A  is  constant  (which  seems  to  be 
true  within  reasonable  error  bounds) ,  the  conversion  to  a  time-dependent 
function  is  accomplished  by  expressing  d  as  a  function  of  time  using  the  vortex 
trajectory  data  of  Fig.  8.  The  next  step  is  the  determination  of  the  residual 
velocity,  Vp,  which  must  be  vectorially  added  to  before  a  comparison  can  be 
made  with  the  probe  data.  He  have  done  this  in  two  different  ways: 

1)  The  first  way  to  determine  involves  very  young  vortices  (about  50* 

old) .  For  these  it  was  assumed  that  Vp  was  given  by  the  vortex  trajectories 
(Fig.  8).  The  fitting  process  always  commenced  when  the  vortex  core  hit  the 
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probe  and  ended  when  the  following  blade  passed  over;  this  assured  the  simplest 
possible  flow  field,  as  there  would  be  vortices  on  only  one  side  of  the  probe 
and  minimal  Influence  of  vortex  sheets  and  blade  bound  vortlclty.  Figure  9  shows 
some  typical  comparisons  of  probe  data  with  the  fitting  expression.  This  figure 
shows  the  vortex  velocity-time  traces  for  pitch  settings  of  8°  and  12**.  It  Is 
seen  here  that  the  fitting  curve  provides  a  good  match  to  the  data  outside  of 
the  Immediate  core  region.  Furthermore,  the  vortex  strength  Is  very  close  to 
the  maximum  blade  bound  circulation. 

2)  A  second  means  to  determine  was  required  In  analyzing  older 

vortices  (about  210**  old).  The  flow  Is  more  complex  In  this  case,  as  the  probe 
always  lies  between  two  vortices  In  the  fitting  region,  and  the  expression  for 
the  vortex-induced  velocity  Is  correspondingly  complicated.  In  fact,  V^,  for 
this  case  was  determined  using  three  vortices  —  one  outboard  of  the  probe  and 
two  Inboard.  Again,  the  data  were  fit  for  the  time  period  between  a  probe- 
vortex  strike  and  the  subsequent  blade  passage.  It  was  found  that  with  Vg 
determined  by  the  vortex  trajectory  data.  It  was  not  possible  to  obtain  a  good 
fit  of  the  classical  vortex  expression  to  the  wake  data.  Instead  we  found  that 
a  better  value  for  Vg  was  found  by  use  of  the  minimum  measured  velocity  between 
two  vortices.  At  this  point,  the  vortex-induced  velocity  Is  small,  but  not 
zero  (due  to  the  differing  Instantaneous  translation  velocities  of  the  three 
vortices) .  The  minimal  Induced  velocity  Is  calculated  (assuming  some  value 
of  A)  and  subtracted  from  the  minimum  measured  Inter-vortex  velocity  to 
obtain  Vg.  This  task  was  rendered  quite  simple  by  the  fact  that  the  radial 
component  of  these  velocities  turns  out  to  be  very  small  (this  was  checked  by 
calculations  and  measurements  with  a  second  probe) .  Since  the  two  methods 
above  do  not  give  the  same  value  for  the  residual  velocity.  It  Is  clear  that 
Vg  Is  not  a  constant  In  this  case.  We  assume,  however,  that  It  changes  suf¬ 
ficiently  slowly  to  render  the  fitting  process  meaningful.  In  fact,  the  results 
thus  obtained  are  consistent  with  the  young  vortex  data.  Figure  10  shows  some 
typical  comparisons  of  the  older  vortex  data  with  the  fitting  expressions..  This 
figure  shows  the  fittings  for  pitch  settings  of  5**,  8®,  and  12*.  It  Is  seen 
that  the  8*  and  12*  cases  show  vortex  strengths  which  match  the  maximum  blade- 
bound  vortlclty  very  well.  At  5*  pitch,  however,  the  strength  Is  seen  to  be 
considerably  smaller. 

It  seems  from  the  above  data,  which  are  taken  at  a  low  rotor  speed,  that 
the  tip  vortex  develops  Its  full  strength  very  early  In  life  (mainly  before  50*) . 
Although  there  Is  a  fair  amount  of  variability  between  vortices.  It  Is  rather 
striking  that  very  many  vortices  closely  approach  a  classical  Ranklne  vortex 
In  appearance.  Furthermore,  the  vortices  (except  for  the  5*  case)  seem  to  con¬ 
tain  all  of  the  blade  circulation.  This  vortex  strength  and  structure  differs 
markedly  from  the  result  obtained  by  Cook  and  probably  reflects  the  considerable 
differences  In  blade  geometries.  As  rotor  tip  speed  Increases  (Fig.  11) ,  how¬ 
ever,  there  appears  to  be  an  Increasing  departure  from  the  Ranklne  vortex 
appearance.  Nevertheless,  the  nondlmenslonal  vortex  strength  seems  imaffected 
by  tip  speed. 

4.  Comparison  of  Theory  and  Experiment 

In  order  to  Integrate  the  present  wake  and  loading  data  Into  a  believable 
whole.  It  Is  necessary  to  be  able  to  reproduce  the  blade  loading  computationally. 
We  have  chosen  to  do  this  using  A.M.I.'s  lifting  surface  code  [6].  This  Is  a 
very  flexible,  compressible,  lifting  surface  code  which  can  handle  either 
prescribed  or  free  wakes. 
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Initial  efforts  to  compute  the  blade  loading  were  done  using  the  Kocurek 
wake  geometry  [5].  The  resulting  computed  thrust  coefficient  was  too  high  by 
about  20%.  The  next  step  was  to  compute  the  loading  using  the  measured  vortex 
locations  and  strength.  Figure  12  shows  a  comparison  of  the  measured  and 
computed  loading  using  the  measured  vortex  parameters  for  a  collective  pitch 
of  8‘  (the  trajectory  Is  given  by  Fig.  8  and  we  choose  A  ■■  1.0).  The  comparison 
Is  now  considerably  Improved  and  the  thrust  coefficient  Is  over-predicted  by 
less  than  5%.  In  view  of  the  previously  mentioned  uncertainties  In  the  vortex 
trajectory  measurements,  these  computations  were  also  performed  with  the  entire 
vortex  trajectory  perturbed  such  that  at  =  180®,  the  axial  and  radial  perturba¬ 
tions  were  ±0.025R.  The  results  derived  from  all  possible  cond>lnatlons  of  these 
axial  and  radial  changes  fill  the  shaded  area  In  Fig.  12.  That  the  above 
measured  and  computed  results  are  roughly  centered  on  this  shaded  region  Indicates 
that  for  this  case  the  measured  trajectories  are  fairly  accurate.  However,  the 
best  comparison  with  the  measured  loading  occurs  when  the  vortex  radial  location 
(at  y  =  180®)  Is  Increased  (that  Is,  the  contraction  Is  decreased)  by  0.025R. 

The  Identical  situation  was  found  to  occur  In  computations  of  the  12®  pitch 
cases;  l.e.,  the  best  comparison  occurred  when  the  radial  vortex  location  was 
Increased  by  0.025R  over  the  measured  value  (Fig.  13).  For  the  5®  collective- 
pitch  case,  the  situation  was  a  little  different  In  that  a  reasonable  comparison 
of  computation  and  loading  data  could  not  be  made  until  the  vortex  strength  was 
reduced  to  A  =  0.75.  In  this  case,  the  vortlclty  which  would  otherwise  have 
been  In  the  tip  vortex  was  now  Included  In  the  vortex  sheet  model.  (For  a 
complete  description  of  the  assumed  vortex  sheet  model  see  Ref.  [6]).  This 
result  Is  consistent  with  the  measured  vortex  strength  and  gives  the  comparison 
shown  In  Fig.  14. 

The  previous  comparisons  have  been  made  at  low  tip  Mach  numbers.  The 
lifting  surface  code  used  should  be  applicable  to  predict  the  spanwlse  and  chord- 
wise  loading  up  to  the  onset  of  supercritical  flow.  Beyond  this  point,  linear 
aerodynamics  are  not  applicable  on  the  blade  and  a  more  complete  flow  description 
Is  required.  As  a  preliminary  evaluation  of  the  high-speed  flow  data,  two- 
dimensional  computations  were  made  of  the  flow  at  the  80%  radial  station.  This 
was  done  using  Holst's  full-potential  code  [17].  In  order  to  perform  this 
computation,  an  angle  of  attack  Is  required.  Since  the  region  of  supersonic 
flow  Is  localized  (l.e.,  limited  to  the  Immediate  vicinity  of  the  upper  blade 
surface) ,  It  should  be  possible  to  find  the  angle  of  attack  using  the  linear 
lifting  surface  code.  Of  course,  the  lifting  surface  code  requires  the  measured 
vortex  location  and  strength  as  mentioned  previously.  With  the  angle  of  attack 
obtained  thereby,  the  Holst  code  produced  the  results  shown  In  Fig.  15.  This 
figure  shows  two  computed  results  —  an  Invlscld  result  and  one  with  a  viscous 
ramp-boundary  layer  model  [18] .  It  Is  seen  that  a  shock-boundary  layer  Inter¬ 
action  model  Is  very  necessary  and  In  this  case  very  effective. 

5.  Concluding  Remarks 

The  present  study  was  Intended  as  a  benchmark  to  aid  In  the  development 
of  hover  performance  codes.  The  goal  was  to  obtain  simultaneous  measurements 
of  blade  load  distribution  and  tip  vortex  geometry  and  strength  using  fairly 
standard  techniques.  In  spite  of  some  uncertainties  (due  mainly  to  wake 
unsteadiness) ,  lifting  surface  computations  show  that  the  load  and  wake  measure¬ 
ments  are  generally  consistent  with  each  other. 

The  main  conclusions  from  this  study  are: 

1.  The  Cook  vortex  measurement  technique  seems  to  be  quite  effective  for 
two-bladed  rotors. 
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2.  At  low  rotor  speeds,  an  untwisted,  untapered,  double-bladed  rotor 
produces  tip  vortices  which  can  closely  resemble  a  classical  Ranklne  vortex. 

Except  for  the  lowest  pitch  settings,  this  vortex  strength  closely  approaches 
the  maximum  blade  bound  circulation.  At  higher  tip  speeds,  the  Inner  vortex 
structure  appears  Increasingly  nonclasslcal;  however,  the  strength  Is  unaltered. 

3.  It  Is  not  possible  to  predict  the  blade-spanwlse-load  distribution 
without  accurate  vortex  location  and  strength  data.  The  present  measured  vortex 
location  data  were  significantly  different  (for  presently  unknown  reasons)  from 
the  classical  data  In  the  literature.  However,  these  measurements  were  Indispen¬ 
sable  to  obtaining  a  reasonable  comparison  of  theory  and  experiment. 

4.  For  the  present  rotor  and  speed  range  tested,  the  onset  of  transonic 
flow  was  found  to  have  no  effect  on  the  spanwlse  loading  distribution  and  the 
vortex  trajectories.  The  chordwlse  loading  Is  profoundly  altered  by  the  transonic 
flow  and  can  only  be  simulated  by  nonlinear  aerodynamic  tectmlques  which  employ 

a  shock-boundary  layer  interaction  model. 
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Fig.  1.  The  model  and  experimental  set-up. 
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Measured  pressure  distributions.  Collective  pitch,  6, 


TEST  DATA,  6g  =  8° 

a  °  1250  rpm.  C-f  •  0.00460 
n  -  2060  rpm,  Cj  =  0.00461 
n  ~  2500  rpm.  Ct  >  0.00464 
REGION  OF  LOADING  VARIATION 
DUE  TO  RPM  CHANGES 


.6  .8 

r/R  (RADIAL  STATION) 


Fig.  5.  Effect  of  rotor  speed  on  blade  span  loading. 
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Fig.  9.  Typical  vortex  velocity-time  trace  and  1/R  curve  fit  for  various 
pitch  settings.  Vortex  age  ■  50”  (nominal).  0  ■  1250  rpm. 


C(  (SECTIONAL  LIFT  COEfFICIENT)  ^  Cg  (SECTIONAL  LIFT  COEFFICIENT) 


Fig.  14.  Comparison  of  measured  and  computed  loading. 


Comparison  of  measured  and  computed  chordwise  pressure  distribution. 
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